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A quantum-mechanical framework is set up to describe the full counting statistics of particles 
flowing between reservoirs in an open system under time-dependent driving. A symmetry relation 
is obtained which is the consequence of microreversibility for the probability of the nonequilibrium 
work and the transfer of particles and energy between the reservoirs. In some appropriate long-time 
limit, the symmetry relation leads to a steady-state quantum fluctuation theorem for the currents 
between the reservoirs. On this basis, relationships are deduced which extend the Onsager-Casimir 
reciprocity relations to the nonlinear response coefficients. 
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S ' T. INTRODUCTION 

Quantum systems can be driven out of equilibrium by time-dependent perturbations, by interaction with reservoirs 
at different chemical potentials or temperatures, or by combining both. In the latter cases, the quantum system is 
open and currents of energy and particles are flowing across the system. Such processes take place in mesoscopic 
electronic conductors as well as in chemical reactions. These nonequilibrium processes can be characterized by the 
relations linking their currents to the differences of chemical potentials. These relations may be linear in the case of 
Ohm's law, but are typically nonlinear, which defines the nonlinear response coefficients. 

Alternatively, the full counting statistics of the particles transferred between the reservoirs can be considered. This 
statistics aims to characterize the transfers of particles in terms of the functions generating all the statistical moments 
£f) • of the fluctuating numbers of particles. The knowledge of this generating function gives access not only to the 
conductance and the noise power but also to higher-order moments and thus to the properties of nonlinear response. 
The full counting statistics has attracted considerable theoretical interest and is also envisaged in experimental 
measurements dUli. After the pioneering work of Ref. [H , several methods have been developed in order to obtain 
the full counting statistics in mesoscopic conductors. One of them is based on Keldysh Green's function formalism, in 
which an expression for the generating function has been obtained within a semiclassical approximation [H, 0, @, Bpl- 
. ,-h , The full counting statistics can also be obtained on the basis of quantum Markovian master equations describing 
the fluctuations of the currents as well as in terms of stochastic path integrals 11 1. The generating function 



obtained in the approaches using the semiclassical approximation or the Markovian master equation has been shown 
to obey a symmetry relation as the consequence of time reversibility In nonequilibrium statistical mechanics, this 
relation is known as the fluctuation theorem which has been established for several classes of systems. These latter 
are either time- independent deterministic [HI, EiL EEJ or stochastic systems sustaining nonequilibrium steady states 
[H, E3, EH EH H3j EH, US, EH, or time-dependent Hamiltonian or stochastic systems, in which case the fluctuation 
theorem is closely related to Jarzynski nonequilibrium work theorem [24], El, Ell . Similar symmetry relations have 
been considered for continuous-time random walks 12711 . Quantum versions of the fluctuation theorem and Jarzynski 



nonequilibrium work theorem have also been obtained [2j| |2g, |3d, lM IM IM ^ 

Moreover, a further relationship has recently been proved for time-dependent quantum Hamiltonian systems 44j, 
allowing the derivation of Green-Kubo formulas and Onsager-Casimir reciprocity relations for the linear response 
coefficients [H [H, 51 51 ■ 

An open question is to bridge the gap separating the time-dependent situations from the nonequilibrium steady 
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states which are expected to be reached in the long-time limit. The problem is to deal with nonequilibrium steady 
states without relying on the semiclassical or Markovian approximations, or on the neglect of the energy or particle 
content of the subsystem coupling the reservoirs. 

In the present paper, our purpose is to prove the fluctuation theorem for the currents in open systems obeying 
Hamiltonian quantum dynamics and sustaining nonequilibrium steady states in the long-time limit. We start by 
considering a time-dependent quantum system in contact with energy and particle reservoirs at different temperatures 
and chemical potentials. The amounts of energy and particles which are exchanged between the initial and final times 
are determined by quantum measurements. This framework is similar to the one considered by Kurchan to obtain 
a fluctuation theorem for quantum systems [28]. Here, this framework is extended by taking the initial states as 
grand-canonical instead of canonical equilibrium states, which allows us to deal with transfers of particles between the 
reservoirs. In this way, we obtain an exact relationship which is the consequence of microreversibility for the probability 
of a certain exchange of energy and particles between the reservoirs during the time-dependent external drive. An 
equivalent symmetry relation is obtained for the generating function of all the fluctuating variables. However, these 
symmetry relations are expressed in terms of the temperatures and chemical potentials of the reservoirs. The problem 
is that we need a symmetry relation in terms of the differences of temperatures and chemical potentials, which define 
the thermodynamic forces (also called the affinities) driving the currents across the system. The importance of this 
point has recently been discussed in the review (40| . 

The central contribution of the present paper is the proof that, in the long-time limit, the aforementioned generating 
function only depends on the differences between the parameters of the reservoirs. This proof is carried out by 
obtaining lower and upper bounds on the generating function in terms of a new generating function which only 
depends on the differences of parameters and further functions which are bounded in the long-time limit. Combining 
this fundamental result with the previously established symmetry relation of the generating function, the fluctuation 
theorem is proved for nonequilibrium steady states in open quantum systems. Thanks to this quantum fluctuation 
theorem, the Onsager-Casimir reciprocity relations and their generalizations to the nonlinear response coefficients can 
be inferred [H|,[4i]. 

The plan of the paper is the following. The protocols for the forward and reversed drives of the open system are 
introduced in Sec. [ill The symmetry relations for the probability and the generating function are proved in Sec. IIIII 
In Sec. IIV1 we obtain the quantum fluctuation theorem for the currents in the steady state reached in the long-time 
limit. In Sec. El the consequences of the fluctuation theorem on the linear and nonlinear response coefficients are 
deduced. The conclusions are drawn in Sec. IVII 



II. OPEN QUANTUM SYSTEM AND TIME-DEPENDENT PROTOCOLS 



A. The total Hamiltonian 



We consider a total quantum system composed of a subsystem in contact with several reservoirs of energy and 
particles. Initially, the reservoirs are decoupled. During a time interval T, the reservoirs are put in contact by some 
time-dependent interaction which has the effect of changing the energy and the particle numbers in each reservoir. 
The total Hamiltonian of the system is thus given by 

r 

H(t; B) = J2 H j for t < > C 1 ) 

3=1 
r 

3=1 
r 

= Y^Hj for T<t, (3) 

3=1 

where Hj denotes the Hamiltonian of the j isolated reservoir before the interaction is switched on. During the time 
interval < t < T, the system is submitted to the time-dependent interaction V(t), which describes the coupling of 
the reservoirs by the subsystem. Beyond the final time t — T, the reservoirs arc decoupled into the Hamiltonians Hj . 
We suppose that the whole system is placed in a magnetic field B. 

The observables of the reservoirs are the Hamiltonian operators Hj and Hj , as well as the numbers of particles of 
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several species Nj a with a = 1, 2, c. The total number of particles of species a is given by 



3=1 



(4) 



where n Q denotes the number of particles of species a in the subsystem between the reservoirs. 

Since the numbers of particles of species a are conserved within each isolated reservoirs before and after their 
coupling, we have that 



j'a\ 







and 



Hj,Nj' a 



0, 



for all j, f — 1, 2, r and a — 1, 2, c. 

We suppose that the Hamiltonian operator H(t; B) has the symmetry 

QH(t;B)e = H(t;-B) 



(5) 



(6) 



under the time-reversal operator 0. This latter is an antilinear operator such that Q 2 — I and which has the effect 
of changing the sign of all odd parameters such as magnetic fields. Equation ([6]) expresses the microreversibility in 
an external magnetic field. The numbers of particles are symmetric under time reversal: 



QN ja Q = N ja and 9n Q 9 



(7) 



In the following, we consider protocols with two quantum measurements at the initial and final times separated by 
a unitary time evolution (see Ref. [40] for a review). 



B. The forward protocol 

The forward time evolution is defined as 

ih^-U F (t; B) = H{t; B)U F (t; B) , (8) 
with the initial condition f7p(G; B) = I. In the Heisenberg representation, the observables evolve according to 

A F (t) =U F (t;B)AU F (t;B) , (9) 
which also concerns the time-dependent Hamiltonian 

H F (t) = U F {t;B)H(t;B)U F {t;B). (10) 
The average of an observable is thus obtained from 

(A F (t)) =trp(0; B)A F (t). (11) 

We note that the dependence on the magnetic field is implicit in these expressions. 

The initial state of the system is taken as the following grand-canonical equilibrium state of the decoupled reservoirs 
at the different inverse temperatures (3j = l/ik^Tj) and chemical potentials fij a : 

p{0 -B) = n ~ g( m = ne-frfo-Sq*"**--^*)] , (12) 

where $>j(B) = — fee?} lnSj(i?) denotes the thermodynamic grand-potential of the j th reservoir in the initial equilib- 
rium state. 

A quantum measurement is performed at the initial time. The system is observed in the eigenstate \^k) of the 
reservoir operators of energy and particle numbers: 

t = 0: Hj\* k ) = e jk \V k ) , (13) 
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At the final time, another quantum measurement is performed and the system is observed in the eigenstate l^;} of 
the reservoir operators of energy and particle numbers: 

t = T: Hj\h) = , (15) 

N ja \*i) = 5,-ai|*i> • (16) 
Accordingly, during the forward time evolution, the following change of energy in the j th reservoir is observed: 

Ae j = Zjl - > (I 7 ) 
while the number of particles of species a in the j th reservoir changes by 

Avj a = Vj a i — Vj a k . (18) 

C. The reversed protocol 

The evolution operator of the reversed process is defined as 

ih^U K (t- B) = H(T - t; B)U R (t; B) , (19) 

with the initial condition ?7r(0; B) — I, and is related to the one of the forward process by the following 

Lemma: The forward and reversed time evolution operators at the final time T are related to each other by 

U R (T;-B) = eul(T;B)Q. (20) 



This lemma is proved by noting that the forward time evolution in the magnetic field B, followed by the operation 
of time reversal, by the reversed time evolution in the magnetic field —B, and finally by time reversal again is equal 
to the identical operator: 

eU R (T;-B)QU F (T;B)=I , (21) 

from which we deduce Eq. (l2"0l) . 

The reversed protocol is supposed to start with the following grand-canonical equilibrium state of the final decoupled 
reservoirs: 

p(T;-B) = 11 ~ = ne-ftfo-S.^.-g.C-B)] > (22) 

3=1 ^' 3=1 

at the same inverse temperatures f3j = \/(k&Tj) and chemical potentials (ij a as in the forward protocol and where 
$j(—B) = — ksTj TnSj(— B) denotes the grand-canonical thermodynamic potential of the j th reservoir in the final 
equilibrium state and reversed magnetic field. 

Similarly to the forward protocol, quantum measurements are performed at the initial and final times to determine 
the changes of energies and particle numbers in the reservoirs. 

III. CONSEQUENCES OF MICRO REVERSIBILITY 

A. The symmetry relation for the probability of the fluctuations 

The probability distribution function to observe the energy (fTT|) and particle transfers (fT5)) during the forward 
protocol is defined as 

p F (Aej,Aiy ja :B) = ^ J| 6 [Aej - (iji - e jk )} ]J 5 [ A ^J" ~ food ~ v j a k)\ 

kl j ja 

x|(* ; (B)|C/ F (T;S)|vI/ fe (i3))| 2 {* k (B)\p(0; B)\* k (B)) . (23) 
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We notice that this function is a probability density because the quantities S(-) are Dirac distributions for both the 
energy and the particle numbers. 

Inserting the expression of the initial density matrix (fT2| . using the Dirac delta distributions to replace the initial 
energies and numbers into the final ones, we find that 

p F (Aej,An a -B) = ^ JJ<J[Ae,- - (iji - e jk )} JJtf [Aj/,- q - (i> ja i - Vjak)] 

kl j ja 

x \($i(B)\U F (T; B)\^ k (B)}\ 2 e" ^ ft ^«^-«*-*,(s)] 

kl j ja 

x\(^i(B)\U F (T;B)\^ k (B))\ 2 e~^^ fa' (-B)] j 

(24) 

where we have introduced the difference of the thermodynamic grand-potential of the j th reservoir as 

A$,= - . (25) 



According to the lemma ([201) , the probability of the transition A: — > Z during the forward process is equal to the 
probability of the transition Z — » fc in the reversed process and magnetic field: 

\(MB)\U V (T;B)\* k {B))\* = |(* z (B)|e£/t (T; -i?)6|* fc (i?)) | 2 

= \(M-B)\U K (T;-B)\* l (-B))\* . (26) 

Substituting this identity in Eq. (|24p and introducing the probability of negative changes in the energies and particle 
numbers during the reversed process as 

PB.(-Aej,-Avj a ;-B) = ^ ±L S [~^ e j ~ ( e jk - tji)] \\S [-Auj a - (uj a k - i>jai)] 

kl j ja 

x\(^k(—B)\Un(T; -B)\^i(-B))\ 2 (&i(-B)\p(T; -B)\^i(-B)) 

(27) 

with the final density matrix (|22p , we obtain the following symmetry relation: 

Pf {Aej , Avj a ; B) = e^ ft( A ^£* w.A^.-A*,) pR (_ Ae „ -Au ja ; -B) . (28) 

If this relation is restricted to the energy change in a single system, this fluctuating quantity is the work W performed 
on the system and we recover the quantum version of Crooks' fluctuation theorem 

p F (W;B) =eP (w - AF) PR (-W:-B) (29) 

with the corresponding difference of free energy AF = F(—B) — F(B) 25]. The relation (j2"8|) extends this result to 
the transfer of particles under the effect of the differences of chemical potentials driving the system out of equilibrium. 

B. The symmetry relation for the generating function 

The generating functions of the statistical moments of the exchanges of energy and particles are defined by 

Gf,r(&, Vj*,B) = J II dAe i dA ^ a e ~ Ej ^ Ae3 " E " p F ,K(At 3 ,Av oa - B) (30) 

ja 

for the forward and reversed processes. The knowledge of these generating functions provides the full counting 
statistics of the process. We notice that the generating function of the forward protocol is alternatively defined as 

G F (^,ri ja ;B) = (e-E i bH&-T; ia Ti i c.Nirt f £ i 6 j ni+T; ia njc.N ia ^ (31) 
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with 

H jF = uI(T;B)H j Ut(T;B), (32) 
N jaF = U F (T;B)N ja U F (T;B) , (33) 

and the generating function of the reversed protocol as 

G R (^,r} ja ;-B) = (e - ^ &fijR-E Ja w-^«.R e Sj «^+E 3 - a \ (34) 



R 

with 

ff iH = UI(T--B)H J U K {T--B) , (35) 
A^ R = E/£(T; -B)N ja U R (T; -B) . (36) 

Taking the Laplace transforms of the symmetry relation (128p , we obtain an equivalent symmetry relation in terms 
of the generating functions: 

<m <,.//„,:/;> = c-M A *> Gr(/3j- — ^j, —f3j[ij a — Tjj a ; —B) , (37) 

in terms of the temperatures and chemical potentials of the reservoirs. As mentioned in the introduction, this 
symmetry relation has not yet the appropriate form because the thermodynamic forces or affinities do not appear. 

IV. QUANTUM FLUCTUATION THEOREM FOR THE CURRENTS 

In this section, we prove that, in the long-time limit, the generating functions entering into the symmetry relation 
(I37p only depend on the differences of the parameters £j and rjj a , leading to the requested symmetry. In the long-time 
limit, a nonequilibrium steady state can be reached between the reservoirs if the coupling remains constant over the 
whole time interval except finite transients. 

A. The theorem 

We consider a situation where two large quantum systems interact through a bounded time-dependent perturbation 
described by V(t). Then, the generator of time evolution of the whole system is given by 

H = H 1 + H 2 + V(t) . (38) 

Hereafter, we assume that V(0) = 0, V(t) = V(T - t) and V(t) = V for t < t < T - t . 

Let Nj a (j — 1,2) be the number of a-particles in the j th large system and assume that [iVj. Q + A^a + n a ,H\ + 
B-2 + Vq] = with a bounded n a (imagine that a quantum dot is located between the two electrodes). 

Since the interaction is symmetric under time reversal V(t) = V(7~ — t), the evolution operator of the forward and 
reversed protocols are identical 

U F (t;B) = U K {t;B) = U(t;B) (39) 
and therefore solutions of one and the same equation: 

i^-U(t;B) = [H 1 +H 2 + V(t)]U(t;B) (40) 
at 

with the initial condition U(0;B) — I and % = 1. For the same reason, the initial and final reservoir Hamiltonians 
are the same, Hj = Hj for all j — 1, 2, r, so that the initial and final density matrices have the same definition 

P(B) = f[e-^-^H a N ja -^(B)} ^ (41) 

where $>j(B) is the thermodynamic grand-potential of the j th reservoir in magnetic field B. Accordingly, the forward 
and reversed generating functions also have the same definition 

G(^,r) ja ;B) = (e~^ C*^-E,„ r, ja N jotF e £^ ,C,flj+£ Ja ni.»JVi.\ ( 42 ) 
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with Eqs. (|32| and (j33)) and where the average (•} is carried out with respect to the density matrix (|4Tj) . 
According to Eq. (|37[) . this generating function has the symmetry 

Gfe t ri ja ;B) = e~^ ^ A *' G(& - ^,-(3 jt i ja - Vja ; -B) , (43) 

in terms of the temperatures and chemical potentials of the reservoirs. 
Our purpose in this section is to prove the 

Proposition. Assume that the limit 

, T} ja ;B) = — lim i In , %C( ; B) (44) 
exists, it is a function only o/£i — £2 and X]\ a — r)2 a : 

Q(€j,Vja:B) = Q(£l Via ~ Via, B) . (45) 



The interpretation of this proposition is that, because of the finiteness of the subsystem and the interaction Vo, 
the energy and particles lost by the left (respectively right) reservoir are transferred to the right (respectively left) 
reservoir within the overwhelming duration to < t < T — to and, as a result, Q becomes a function Q depending only 
on the differences £1 — £2 and r\\ a — rj2 a - We remark that the explicit form of the generating function 

<3(6 - £2, Via - t?2«; B) is given by Eq. (Sty. 

The above proposition implies that 

mi-^Via-Vi^B) = Q[<:,. n ,,:l>>> (46) 

= Q{Pi-^-PiNa-ma\-B) (47) 

= Q(Pl - $2 - 6l + 6) -PlPla + P2^2a ~ Via + V2a\ ~B) (48) 

= Q(A - £1 + £2, -4a ~ Jfto, + V2a\ -B) , (49) 

where we have introduced the affinities: 



A = 0x - j3 2 , (50) 
A Q EE /3 2 ^2a ~ PlUXa 1 for ot = 1, 2, c, (51) 

driving respectively the heat current and the a-particle currents from reservoir 2 to reservoir 1. The result (|49[) is 
obtained by using the definition (|45]) at the line (|46|) . the symmetry (|43|) and the independency of the quantities A$j 
on the time interval T at the line (|4"T|) . again the definition (|43|) at the line (|4"8")l . and finally the definitions of the 
affinities (|50"|) and (|5ip. Hence, we have the 

Fluctuation theorem. T"/ie generating function of the independent currents satisfies the symmetry 

Q({,r) a ;B) = Q{A -t,A a -ri a ;-B) . (52) 



In the particular case where the two systems have the same temperature, (3\ =02, the generating function has the 
symmetry: 

Q{Si,r la -B)=Q{-tA a ~T la -~B) , (53) 

and we recover the symmetry 

Q(0,r] a ;B) = Q(0,A a -r] a ;-B) (54) 

of the generating function of the independent particle currents, which has already been proved elsewhere for stochastic 
processes [23l ]. 

We notice that the fluctuation theorem (|5"2"1) which is here proved thanks to the proposition (|4"4"1) - (|4"5")) reduces to 
the steady-state fluctuation theorem presented as Eq. (104) in the review [40j for vanishing magnetic field, B = 0. 
Accordingly, the proposition (|4"4"|) - (|4l)|) also provides a rigorous proof of such steady-state fluctuation theorems. 
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B. Setting 

In order to demonstrate the above proposition, the time evolution is decomposed into different pieces corresponding 
to the short initial transient over < t < t , the long steady interaction over i < t < T — t , and the final short 
transient over T — to < t < T. We introduce the lapse of time of the steady interaction 

r = T-2t . (55) 

In addition to U(t; B) defined by Eq. {§]), we introduce U\{t; B) as the solution of 

i^U 1 (t;B) = [H 1 + H 2 + V(t -t)}U 1 (t;B) with U 1 (0;B) = I. (56) 

It is then easy to show 

U(t + 2t ; B) = Uxito; B)e- lST U(t Q ; B) = U f U T U, , (57) 

where H = H 1 + H 2 + V Q , U f = Ui(t ; B), U T = &- lRr and U t = U(t ; B). We further note that 

e iH {r+t ) ui(r + 2to] B)e iHoto = e iHoT T f e- iHoT T T e- iHoto Tie iHoto , (58) 

where H = H 1 + H 2 , T f = e lH ° to U f: T T = e iH ° T U T , and T z = e %H ° t °U l . 
Therefore, with the aid of [Ho, Hj] = [Ho, Nj a ] — [Hq, p] — 0, we have 

G(£j,rjj a ;B) 

= (U(t + 2t ; B)t e -^i - 11 > : • n- N i« U(t + 2t ; B) e^j *i B i+T. ia n ia N ja ^ 
= (C/tyt^te-E^i^'-EiaWciV^ UfUrU . eE^fij+EicW-Wia) 

= (T i (-to^TlT f (T)h-^i^-^ r >i aN ^ Tf(T)T r Ti(-to) e^i ei H i+^^^ N i«) , (59) 

where T A (r) = e lHoT T x e-' HoT (A = i or /). 
For later purpose, we introduce 

E & H i + E VjaNja = 2C + 2D, (60) 

3 jot 

2C=^±^H Q + ^ V -^^N Qa7 (61) 

a 

2D = (6 - 6)Aff + ^(ry lQ - r) 2a )&N 0a , (62) 



2 ^ = E^ (^-E^^) > 

i \ « / 



(63) 



where H = H 1 + H 2 , N 0a = N la + N 2a , AH = {H x - H 2 )/2, and AN 0a = (N la - N 2a )/2. Since [C, p] = [D, p] = 0, 
we have 

G(^,r, ja ;B) = (e^I^-toMr f (r)^- 2C - 2D T f (T)T T T t (-t a ) c c+D ) . (64) 

This is our starting point. Note that C and D are Hermitian for real £j and T]j a and that -D is the function only of 
£l - 6 and 7yi a - r] 2a . 



C. Some Inequalities 

Here, for the sake of self-containedness, well-known equalities and inequalities [5(| necessary for the following proof 
are summarized. For an operator X, the operator norm \\X\\ is defined by 



\xy , S up Jm^m (65) 
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and it satisfies: 

Equality 1: For any unitary U, \\W XU\\ = \\X\\. 
Indeed, we find 

_ wxxm _ m^xm _ sup <m> _ 

(<P\<p) w^o (il>\UW\i/>) Wjt0 (ip\ip) 

where we have set \ip) — U\ip). 
Inequality 1: (XWYX) < \\Y\\ 2 (X^X). 

Let {(fa} be a complete orthonormal basis of eigenvectors of p: p\<p a ) — p a \ip a ). Then, because of (ip\X^X\ip) < 

\\x\\ 2 {<p\<p), 

{xWyx) = Y,pM x]y]yx \<p°) 
< J2pAy\\ 2 (^x\^) 

<7 

= \\YfY,P^a\X^X\ Va ) 

a 

= \\Y\\ 2 (X^X) ■ (67) 

Inequality 2: (X^YX) < \\e- A Xh A \\ 2 (Y^Y) where 2A = £\ pj(Hj - *Z a ii ja N ja ). 

Thanks to the cyclicity of the trace, we have the Kubo-Martin-Schwinger (KMS) condition (XY) = 
(e A Ye~ A e~ A Xe A ) for canonical averages (XY) = =Ty(c~ 2A XY) with S = Tre~ 2A . The KMS condition and In- 
equality 1 imply 

{XWYX) = (c A YXc- A c- A X^c A ) = (c A Yc- A c A Xc- A c- A Xh A c- A Y^c A ) 

< \\e- A Xh A f(e A Ye- A e- A Y^e A ) = \\e- A X*e A \\ 2 {y*Y) . (68) 



D. Proof 

Step 1: 

Let Xi = c- c - D T T ri(-t ) c c+D . Then Inequality 1 leads to 

G(^, Vja ;B) = (Xl[e- G -^ J (r)e G + D ]^[e^- D r f {r)e c +^X 1 ) 

< \\e- c - D r f (r) e^fiX^) . (69) 

Since r/(r)^F/(r) = 1, we have 

(Xlx,) = (xlc c+D c- c - D e- c - D e c+D X 1 ) 

= (X 1 t c ^F / (r)tF / (r)e- c ^c- c ^F / (r)tF / (r)e^X 1 ) 

= {Xte c + D T f (T)U- c - D [e c+D T f (T)e- c - D ] [e~ c - D r f (r)U c + D ] e~ c - D T f {r)e c+D X x ) 

= (x 1 t e c+Z3 r / (r)t e - c - D y 1 t y 1 c- c - D r f ( T )e c+D x 1 ) 

< \\Y 1 \\ 2 (X\e c+D T f {T)U- 2C - 2D T f {T)e c+D X l ) 

< \\c- c - D T f (Ty C c + D \\ 2 (xlc c + D T f (T)h- 2C - 2D r f (r) C c + D X 1 ) 

= \\e- c - D T f {T^e c+D \\ 2 G^,ri ja ;B) , (70) 

where Y\ = c~ c ~ D Y /(r)'' c c+d and Inequality 1 has been used. 
Since T f ( T ) = e iH ° T T f e- iH ° r and [H , C] = [H , D] = 0, 

C-D rJ ^a+Du _ \\ a -C-DAH T V , a -iH T JD+D\\ 



||e-°-"I7(T) g u+u \\ = \\e- u - u e* M ° T r f e 



= \\e lHoT e- c - D T fG c+D G -^ T \\ 

= ||e- c - Z3 r / e c+Z5 || , (71) 
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where we have used Equality 1 for the norm. Similarly, ||e~ c ~- c T / (T) t e c ' + - D || = \\e^ c ^ D T^ f e c+D \ 
In short, in terms of 



we have 



Gi&,»&a;B) = (xt x,) = (e c+Z3 r 4 (-to) t rt e - 2C '- 2B r T r 4 (-t ) C c+D ) , 

Gl(tj,Vja,B) 



<G(&,»fe a ;B) < ||e" 



C-D r n+D\\2 

i f e 



Step 2: 

In terms of X 2 = e~ c ~ D ri(— t )c c+D , one has from Inequality 2 

g^^-b) = (e c +^r J (-t )trtc- 2C - 2C r T r J (-t ) e c + c > 

= (c c+I5 r 4 (-to) t c- c - I? c c + I? rte- 2C '- 2I? r T c c + I3 c- c - B r l (-to) c c+I5 ) 
= (4 [e-^r^+Y [e- c - D r T c c+D ] X 2 ) 

< ||c- 4 X 2 t e' 4 || 2 ([c- c '- I3 r T e c + D ] t [e- c ~ D T T e c + D ]) 

< \\c- A c c + D r l (-t y c - c - D c A \\ 2 (e c + D rlc- 2C - 2D r T c c + D ) . 

Conversely, in terms of Y 2 = e^ c ^ D Ti(— to)^e c+D , Inequality 2 leads to 

(c c+D rlc- 2C - 2D r T c c+D ) 

= (c c + B r l (-t )e- c -^e c+I3 r l (-t )trt c - 2 ^- 2 ^r T r l (-t )c c + B c- c - c r 4 (-t )te c + c ) 
= (r 2 t c c + B r 4 (-to) t rte- 2C - 2B r r r 4 (-to)o c+I5 r 2 ) 

< ||c^r 2 t c^|| 2 (e c + c r t (-to) t rtc- 2C - 2B r r r t (-to) c c+D ) 

< ||c^c c + B r l (-t )e- c -^e A || 2 (c c + B r 4 (-t )trt c - 2 ^ 2 ^r T r 4 (-t ) c c+D ) . 

In short, let G 2 (^,r] jce ; B) = {e c + D T\e- 2C - 2D T T e c+D ), then 



\e~ A e c + D r i (-t )c- c ^ D e A \\ 2 



< Gi(£j,rij a ;B) < \\e~ A e c+D T i (—to)^e~ C ~ D e A \\ 2 G 2 (^j,r]j 0: ; B) 



Step 3: 

We set 



2C = — + ^ 2 H + ^2 Vla + 



2 ^2 

where H = H x + if 2 + V and 7V Q = AT la + jV 2a + n Q . Then, in terms of f T = U T e lH ° r , we have 

G 2 (^j, T]j a ; b) = ( c c+D rlc- 2C - 2D r T c c+D ) 

= (c C+D e iHoT T^ c - iH o-r c -2C-2D e iHoTY^ e -iH T e C+D^ 
= (e iH 0T e D e Cft Te -C e -D e -iH T LiH T e D e C e -C e -D e -iH T- 



= (xl 



e iH T e -D e -C e C e D e -iH T 



e iHoT e -D e -C e C e D e -iH T 



e lH ° T e- D e- c f T e c e D e- lHoT ) 
t r - 



e iH T e -D e -C e C e D e -iHoT 



x 3 ) 



< \\^ G - D c - c c c <P G - lii ^f(x\x z ) 



>^^^\\\X\X 3 ) 



(72) 



(73) 



(74) 



(75) 



(76) 



(77) 



(78) 



where X 3 = e iH ° T e- D e- c t T e c e D e- iH ° T and [H ,C] = [H Q ,D] = 0, Inequality 1 and Equality 1 have been used. 
Because of e~ c f T e c = e~ lHoT T T c- c c c c lHoT , one has X 3 = e- D T r e- c e c e D and, thus, 



G 2 (£ j ,r ]ja ;B) < \\c- D C - c c C e D \\ 2 (e D e c e' C Tle- 2D T T e- C e c e D ) 



(79) 
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Furthermore, Inequality 2 gives 

(e D e c e- C Tte- 2D T T c- C 'e c c D ) = ( 



c -^ e - c e c e D 



, D rU- 2D r T e D 



c -D e - c e c e D 



c -D e -c e c e D 



c- D e-°e c e D 



<\\e~ 

Thus, G 3 {€j,Vja;B) = (e D Tle- 2D T T e D ) satisfies 

G 2 (^ Vja ;B) < \\e- D e- c e G e D \\ 2 \\e- A 
Conversely, we have 

G3&,Vja;B) = ([e^e-^e^^^e^-^rte-^r.e 

t 



e A \\ 2 (e D Tle- 2D T T e D ) . 



e A \\ 2 G^ 3 , ma] B) . 



(80) 



(81) 



- c e c e D 



c -D e -c e c e D 



< e" 



c -D e -C e C e D 



c -D e -C e C e D 



e A \\ 2 {e D e c e- c Tle- 2D r T e~ c e c C D ) 



t 



C A || 2 <4*3> 



(82) 



where we have used Inequality 2. Let Y3 = c lH ° T c D c c e c e D e lH " T , then Inequality 1 and Equality 1 lead to 

it 



(4*3) = <*] 



e iH T e -D e -C e C e D e -iH T 



Y 3 % 



e iH T e -D e -C e C e D e -iH T 



< \\Ys\\ 2 (Xl 



e iH T e -D e -C e C e D e -iH T 



e iH r e -D e -C e C e D e -iH T 



x 3 ) 
x 3 ) 



G 2 {£j,ilj a ;B) 



Thus, 



= \\c lH ° T c^ D c~ c e c e D e~ iH ° T ^ 2 
= \\c- D c- C e c e D \\ 2 G 2 (^,r l3a ;B) . 

G 3 (£,j,Vj a ;B) 

\\c- A [ c -- D c- c e c e i3 ] t e A || 2 ||e- D e- <5 e c e I) || 2 



< G 2 (£j,Vja,B) . 



(83) 



(84) 



Step 4: 

From Steps 1 to 3, in terms of 
L = 



\c- c - D T)c c + D \\ 2 \\e- A e c + D T l {-t )e- c - D e A \\ 2 \\e- A [e~ D e- c eCe D ] 1 e A \\ 2 \\c- D e- e e c e D \\ 2 ' 



K = ||e- z, e- < W , || 2 ||e- A 



c- D e- c e c e fl 



e^|| 2 ||e- A e c + D r 4 (^ )te- c - I3 e A || 2 ||c- c '- I5 r / c c+D \\ 2 , 



we have 



L G 3 (£j, T]j a ; B) < G(£j,r)j a ; B) < K G 3 (^,T, ja ; B) . 



(85) 



(86) 



(87) 



Note that the constants L and K are independent of r and that G 3 (£j , r]j a ; B) is a function only of £1 — £2 and 
Via — V2a since the operator D depends only on them. 



Step 5: 
If 



<3(£j , Vj a \B) = - lim — In Gfa , r/ ja ; B) 



(88) 
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exists, one has 



Q(€j,Vj<x]B) = - lim ^hxG(^j,Tf ja ;B)+ lim ^lnL 

< - lim — In G 3 , ry Q ; S) 

< - lim i In , %a 5 B) + lim ^ In ^ = , »fca ; S) ■ (89) 

In short, we have shown: 

Qi&VjalB) = - lim ilnG 3 &,»fe Q ;B) = - lim - m^lVV^IVe ) . (90) 
The left-most term only contains D, which is a function only of £i — £2 and rj\ a — 772a, or 

Q(^,Vj a ;B) = Q(£t -£ 2 ,77i Q ~T] 2a ;B) . (91) 

We would like to remark that, even when the system has very long but finite recurrent times, the quantities Q and Q 
can be introduced and the Proposition is expected to hold with errors of order of 1/T. Firstly, in such a case, the ratio 
— In G/T would converge to a definite value Q provided T is sufficiently longer than the relaxation time but shorter 
than the recurrent time. The same would be valid for Q. Secondly, even in such a case, the quantities appearing 
in l|85p and ([55]) are bounded by constants independent of the reservoir volumes fi if the interaction V(t) and the 

subsystem particle numbers n a have finite norms. Indeed, to) = e" ,J?ot °T exp —i f*° e lH ° s V(s)e~ lH ° s ds e lHot ° 

where 'T exp' is the time-ordered exponential, and its norm is bounded by exp [to sup t || V(t)||] [50| . Then, the product 
e c+D Ti(— to)e~ G ~ D is Sl-independent and has a finite norm because e~ c ~ D is 0(6^^) if c c+D is 0(e ±f2 ). By similar 
arguments, all the norms in Eqs. ([55)) and (fBT)]) are found to be ^-independent and, thus, the difference Q — Q is of 
the order of 1/T instead of fi/T. Accordingly, the equality (|9"Tj) is obtained in the limit T — ■> 00. 



E. Generalization 



The previous results can be generalized to the case of r > 2 reservoirs. In this case, the proposition (|45|) is that the 
generating function is a function 

Q(Z j ,r ]ja ;B) = Q(£ j ,f) ja -,B) , (92) 
depending only on the independent parameters: 

k=l 
1 r 

fjja = Vja ~ -Z^Vka , (94) 
k=l 

with j = 1,2, ...,r — 1. The proof is similar as in the case r = 2 with the operators: 

fc=l V a J 



fc=l 



replacing Eqs. ([BT}-((B2]), where i? = X)*=i and ^Oq = Z)/t=i Nka - 
In the general case, the fluctuation theorem should read 



(95) 



Q(£j , Vja \B)=Q (A 3 o - ij , Aj a - fjj a ; —B) , (97) 
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in terms of the independent affinities 

A j0 = &-£X>, (98) 
fc=i 

1 r 

Aj a = -f3j(J,j a + -^PkHka , for a = 1, 2, c, (99) 

with j = 1, 2, r — 1. 

V. SYMMETRY RELATIONS FOR THE RESPONSE COEFFICIENTS 
A. Fluctuation theorem and response coefficients 

If we gather the independent parameters and affinities in the case of r = 2 reservoirs as 

a = {&-6,m a -??2cj , (ioo) 

A = {A ,A a } , (101) 

or in the general case of r > 2 reservoirs as 

A = (102) 

,4 = {A j0 ,A ja } , (103) 

with a = 1,2, c and j = 1, 2, r — 1, the fluctuation theorem (|52|) reads 

Q(A,vl;B) =Q(A-A,vl;--B) , (104) 

where we have explicitly written the dependence of the generating function on the affinities defining the nonequilibrium 
steady state. 

The idea is to differentiate successively the fluctuation theorem with respect to both A and A to obtain symmetry 
relations for the linear and nonlinear response coefficients as well as further coefficients characterizing the statistics 
of the current fluctuations [2l[ . 

On the one hand, the mean currents can be obtained from the generating function and, on the other hand, expanded 
in powers of the affinities: 

J a {B) = ^Q-(Q,A; B)=J2 L a AB)A p + \J2 M a M B ) A P A ~< + lYl N a ,^ s (B)ApA y A B + • ■ • , (105) 

P /3,7 /3,7,<5 



which defines the response coefficients: 



L a , /3(B) = ^(0,0;*), (106) 
Ma,/3-f(B) . dXaB £ BAr mB), (107) 

^<*> s d x ad IXw °^ Bh (io8) 



around the state of thermodynamic equilibrium. 

We notice that, if we set A = in the fluctuation theorem (l52l). we obtain the identities 



Q(0,A;B) = 0, (109) 
Q(A,A:-B) = 0. (110) 
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The former is a condition of normalization and the latter a condition of global detailed balancing which is a consequence 
of the fluctuation theorem ([52]) but may be assumed for itself as a weaker property than the fluctuation theorem [5l[ . 
On the other hand, the generating function of the cumulants of the fluctuating currents at equilibrium satisfies 

Q(X,0;B) = Q(-X,0;-B) , (111) 

obtained from the fluctuation theorem (|52p at the equilibrium A = 0. 

We start by differentiating the fluctuation theorem with respect to the generating parameters {X a } and also the 
affinities {A a } to get 

^L{\,A;B) = -?Q-(A-\,A;-B), (112) 
{\A-B) = ^-(A-X,A--B) + ^-(A-X,A;-B). (113) 



dA a ' dX a ' dA a 

We set A = in Eq. (fTT3j) . use the conditions (fT09|) - pT0)l . and set A = to get 

dQ 



(0,0; B) = 0, (114) 

dX a 

<>Q (0,0; B) = 0, (115) 



dA a 

which shows in particular that the mean currents vanish at equilibrium. 

B. Onsager-Casimir reciprocity relations and Green-Kubo formulas 

Now, we differentiate Eq. (|112[) with respect to to obtain 

d 2 Q d 2 d 

- (\,A;B) = ^-(A-X,A;-B). (116) 



dX a dXj3 ' ' dX a dXp 
Setting A = A = 0, we find the identity 



d 2 2 O 

" -(0,0 ;J B) = -f^-(0,0;-S), (117) 



dXadXp dXadXp 

for the second-order cumulant of the current fluctuations at equilibrium. 
If we differentiate Eq. (|112[) with respect to A@, we get 

d 2 Q d 2 Q d 2 Q 

W {X,A; B) = — ° ^ (A - X,A; -B) - ° % (A - X,A; -B) , (118) 



dX a dAfj dX a dXp dX a dAp 

which reduces to 



L a , (B) = - JZ (0, 0; B) - L a , p {-B) , (119) 

OA a OAp 



for A = A = 0. We recover the formulas of Green-Kubo type in the case a = j3: 



L a , a (B) = -\^(0,0;B). (120) 



The differentiation of Eq. (|113p with respect to Ap leads to 



d 2 d 2 d 2 

( - -(X,A;B) = -^- {A -X,A;-B) + -^-(A-X,A;-B) 



dA a dAp dX a dX/3 dX a dAf, 



2 O d 2 

' -(A-X,A;-B)+ ° f A (A-X,A;-B). (121) 



dA a dXf} y ' ' ' dA a dA, 
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Using Eqs. (fT09|) and (fTTOf in the limit A = A = 0, we have 

= if 2 ?x + L <*.f}{-B) + Lp, a {-B) . (122) 

Combining with Eq. (|119[) , we finally find the Onsager-Casimir reciprocity relations: 

L a , p (B) = L p , a {-B) . (123) 

We notice that Eq. (|119[) leading to the Onsager-Casimir relation requires the link established by the fluctuation 
theorem (f52|) between the variables A and A and does not result from Eqs. (|109[) . (| 1 1 Of) . and (jllip alone. 

C. Symmetry relations at 2nd order 

We proceed in a similar way to obtain relations at second order. The differentiation of Eq. (|116|) with respect to 
A 7 gives 

—^—(X,A;B) = -—^—(A-X,A;-B), (124) 

OXaOApOX-y OXadXpOX^ 

which reduces to 



' r ' Q -(0,0;i?) = - (0,0; -B), (125) 



dX a dXpdXj ' dXudXpdXy 
for A = A = 0. 

On the other hand, the differentiation of Eq. with respect to A y gives 

V (A, A; S) = ax - A,i4; -B) + fl , (ii - A, A; -B) . (126) 



dX a dXfjdAj ' ' dXadXfjdX-y ' ' dX^dX^dA^ 

Here, we introduce the coefficients 

R ^ B ^-wS^ {oABh (i27) 

which characterizes the nonequilibrium response of the second cumulants of the current fluctuations. Setting A = A = 
in Eq. (|126|) leads to the relation 

R aM {B)-R a ^(-B) = dX ® d ® dx (0,0; B) . (128) 
Now, if we differentiate Eq. (1118|) with respect to A 1 , we get 

-(\,A;B) = - — d 3< ? ( A _ A>it ._ g )_ (A- A, A; -f?) 



dXadApdA^ dX a dXpdX~f dXadXpdA^ 

<rQ (A-X,A; -B) - — gg- (A -A, A; --B) 



9A Q 9 J 4 ( g9A 7 ' dXadAfjdA^ 



(129) 



whereupon we find for A = A = that 



M a ,^(B) + M at ^(-B) = R af}>1 {-B) + R ait0 (-B) + m % Q (0,0; fl) , (130) 

OAo.G'AjgG'A.y 



involving the second-order response coefficient (|107|) . 
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We end with the differentiation of Eq. (|121[) with respect to A 1 to obtain 

° U} -(A,A;fl) = -4^(A-X,A;-B) + —4^ r (A-X i A;-B) 



dAadApdA^ ' dXadXpdX^ " dX a dXpdA^ 



' TQ -(A X,A; —B) + - (A - X,A; -B) 



dXadApdXy ' dXadApdA^ 

oAaOApoX^ oAaOXfjoA^ 

d3Q .(A-X,A;-B) + —^——(A-X,A;-B) 



dA a dA p dX 1 y ' ' ' dAadApdAy 



(131) 



Taking A = j4 = therein, we deduce 



M Q , / 3 7 (B)+M /3l7a (S)+M 7) ^(S) = i2 a/ , i7 (B)+JJp 7 , a (B)+il ra|/ 3(B) - - (0,0; g) . (132) 

We notice that this relation is the consequence of the weaker condition of global detailed balancing (|110[) alone and 
could hold even though the fluctuation theorem (p)2"|) does not as recently shown in Ref. [5l| where the versions of 
Eq. (|132p . which are (anti)symmetrized with respect to the magnetic field, appears with the notations M a ^ — 

(ksTfG^, R a0 , 7 = ksTS^ n , and gx Jgg ^ (0, 0; B) = C^. We point out that Eqs. and G3D) are 

on the other hand consequences of microreversibility and the stronger fluctuation theorem, as it is the case for the 
Onsager-Casimir reciprocity relations. 

Adding Eq. (|132[) to the same equation with — B instead of B and using Eqs. (|130p . we moreover infer that 

<! " <} (0,0;B)=Q, (133) 



dX a dXfjdX-y 

whereupon we finally find the symmetry relations: 

R a /3,j(B) — R a i3 ri ( — B) , (134) 
M a , 01 {B) + M aifh (-B) = Rap„(B) + R ai , f j(B) , (135) 
M atfh (B) + M p „ a (B) + M 1>a0 {B) = Ra P „(B) + R 01 , a { B ) + &raA B ) ■ ( 136 ) 

If the magnetic field vanishes B = 0, Eqs. (|135[) reduces to the symmetry relations 

M ai/37 (0) = i [R a p n (0) + R ai ,p(0)] , (137) 
which were previously deduced as consequences of the fluctuation theorem [H, 5§] . 

D. Symmetry relations at 3rd order 

Beside the third-order response coefficient ()108|) . we introduce the coefficients with 

S ^ B ^ 8xJXdA s ^ ^ 



and 

T^s(B)^- dX J^ dAg (0AB), (139) 
characterizing the nonequilibrium responses of respectively the second and third cumulants of the current fluctuations. 
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We continue the deduction by further differentiating the relations of the previous subsection with respect to or 
A s at A = A = 0. We find successively from Eq. (|124j) : 

UiQ -(0,0;fl)= fl , Si 5, (0,0; -B) (140) 



and 

fromEq. |[T25|) : 
from Eq. 1(129)) : 



d\ a dXpdX 1 d\s ' d\ a d\/3d\~ ( d\$ 



W*)+ W~*) = - dx JXdXs { ° AB) ' (W1) 



T a p,ys{B) — T a f3^$(—B) — S a 0-y,s(B) — S a ps,f(—B) , (142) 

N a ,p lS {B) + N a ,/3 lS (-B) = T a p <l5 (-B) + T aiJjS {-B) + T aS . Pl {-B) 

— S a /3f,s{—B) — S a 0S,f(—B) — S aj s,p{—B) 
d 4 Q 



dXadXpdX^dXs 



(0,0; B), (143) 



and from Eq. (pl|) : 



N a ,(3 lS {B) + N 0ty5a (B) + N lM(3 {B) + Ns, a01 {B) 

T a p,ys(B) + T aj ^s(B) + T a 5fiy(B) + Tp 1:a $(B) + Tp$^ ai (B) + T 7 s, a /3(B) 

S a 0y, S(B) - SpyS,a{B) - S 7 Sa,p(B) - Ss a p n {B) 

d 4 Q 



dXadXpdXydXs 



(0,0; B). (144) 



This last relation is the consequence of the weaker condition of global detailed balancing in the same way as Eq. 
In the case of a vanishing magnetic field B = 0, Eqs. (|14ip and Q143P reduce respectively to 

^°» = -l ww P;0) (145) 

and 

N a ,fhs(0) = \ [1^,^(0) + T ayJ3S (0) + T aSj0y (O) - S a07jS (O)] , (146) 
which has been obtained elsewhere as a consequence of the fluctuation theorem [49| . 



VI. CONCLUSIONS 



In this paper, a fluctuation theorem for the currents has been proved for open quantum systems reaching a nonequi- 
librium steady state in the long-time limit [52|. In the considered protocol, the heat and particle currents are defined 
in terms of the exchanges of energy and particles between reservoirs, as measured at the initial and final times when 
the reservoirs are decoupled. 

We start from a general symmetry relation for the generating function of the exchanges of energy and particles, 
which is the consequence of the microreversibility guaranteed by the measurement protocol. This symmetry relation 
is expressed in terms of the temperatures and chemical potentials of the reservoirs. However, the fluctuation theorem 
for the currents requires a symmetry with respect to the thermodynamic forces or affinities which are given in terms 
of the differences of temperatures and chemical potentials. 

We show that this symmetry indeed holds by proving that, in the long-time limit, the generating function only 
depends on the differences between the parameters corresponding to the different reservoirs. Accordingly, the gener- 
ating function of the currents has the symmetry of the fluctuation theorem with respect to the affinities. A rigorous 
proof is thus established for such steady-state fluctuation theorems as considered in the review [40l |. 
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As a consequence, the Onsager-Casimir reciprocity relations can be obtained for the linear response coefficients 
from the fluctuation theorem. Furthermore, generalizations of the reciprocity relations to the nonlinear response 
coefficients can also be deduced. 
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